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Abstract 

Theoretical investigations of surface-state electron dynamics in noble met- 
als are reported. The dynamically screened interaction is computed, within 
many-body theory, by going beyond a free-electron description of the metal 
surface. Calculations of the inelastic linewidth of Shockley surface-state elec- 
trons and holes in these materials are also presented. While the linewidth 
of excited holes at the surface-state band edge (kii = 0) is dominated by a 
two-dimensional decay channel, within the surface-state band itself, our calcu- 
lations indicate that major contributions to the electron-electron interaction 
of surface-state electrons above the Fermi level come from the underlying bulk 
electrons. 
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1. INTRODUCTION 



Shockley surface states are known to exist near the Fermi level in the T-L projected bulk 
band gap of the (111) surfaces of the noble metals Cu, Ag, an Au [0-0]. The wave functions 
of these crystal-induced surface states 0,|5| are localized near the surface and decay into the 
solid, in contrast to the Bloch waves propagating into the bulk. Hence, these states form 
a quasi two-dimensional (2D) electron gas, which overlaps in energy and space with the 
three-dimensional (3D) substrate, and represent a playground for lifetime investigations. 

Various techniques have been used to measure the lifetime broadening of surface-state 
electrons and holes near the Fermi level. High-resolution angle-resolved photoemission spec- 
troscopy (ARPES) has been used to investigate the lifetime of surface-state holes |^-|TT|. In 
contrast to this technique, which is restricted to lifetime measurements of occupied states, 
it has been demonstrated that scanning tunneling spectroscopy (STS) offers the possibility 
to measure the lifetime of long-living surface states above and below the Fermi level fl2|-|14"|| . 
The lifetime of excited holes at the band edge of the partially occupied Shockley surface 
states on Ag(lll), Cu(lll), and Au(lll) was accurately determined by Li et al [0 and by 
Kliewer et al |L3| with the use of STS. Recently, STS has been used to measure lifetimes of 



hot surface-state and surface-resonance electrons, as a function of excess energy ||14j| . 

On the (111) surfaces of Cu, Ag, and Au, the Shockley surface state at the center of the 
surface Brillouin zone (kii = 0) lies just below the Fermi level, with E—Ep = —445, —67, and 
— 505 meV, respectively E and Ep representing the surface-state energy and the Fermi 
level. The dispersion relation of these states is displayed in Fig. 1, which shows that they 
are free-electron-like with effective masses of 0.42, 0.44, and 0.28 m e (m e is the free-electron 
mass) that account for the potential variation parallel to the surface. This figure clearly 
shows that the decay of Shockley surface-state electrons and holes in these materials may 
proceed either through the coupling with bulk states (3D channel) or through the coupling, 
within the surface-state band itself, with surface states of different wave vector parallel to 
the surface (2D channel). 
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While 3D free-electron calculations of the decay rate of excited holes at the surface-state 
band edge (kn = 0) of the noble metals are known to predict significantly longer lifetimes 
than those observed with either ARPES [|6|-|TT|| or STS [|1^,|T^], recent measurements |nj 
have shown that, in the case of surface-state electrons above the Fermi level, the exper- 
imental values are comparable to the calculated lifetimes of bulk free-electrons with the 
same energy. That the decay of surface-state holes is dominated by 2D electron-electron 
interactions, but screened by the underlying 3D electron system, has been demonstrated 
recently |T3|JT5|] , showing an excellent agreement with experiment. In this paper, we extend 
these theoretical investigations to the case of hot surface-state and surface-resonance elec- 
trons in Cu(lll), and show that major contributions to the electron-electron interaction of 
surface-state electrons above the Fermi level come from the underlying bulk electrons. We 
also present self-consistent calculations of the screened interaction, and investigate the role 
that both the partially occupied 2D surface-state band and the underlying 3D electron gas 
play in the relaxation mechanism. 

The contribution to measured lifetimes of surface-state electrons in the noble metals 
arising from electron-phonon interactions has been discussed in |T3|]. Here we focus on the 
investigation of the energy-dependent inelastic lifetimes that are due to electron-electron 
interactions. 



2. THEORY 

We take an arbitrary Fermi system of interacting electrons and consider an excited 
electron or hole interacting with the Fermi sea. The electron-hole decay rate, i.e., the 
probability per unit time for the probe particle to scatter from an initial state ipi(r) of 
energy E± to some available final state ipf(r) of energy Ef (\Ef — E F \ < \Ei — E F \), by 
carrying the Fermi gas from the many-particle ground state to some excited many-particle 
state, may be obtained by using the golden rule of time-dependent perturbation theory ||16|| . 
By keeping terms of first-order in the screened interaction, one finds (we use atomic units 
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throughout, i.e., e 2 = h = m e = 1): 

r- 1 = 2^/dr/dr>*(r)^(r') 

/ 

xlm [-W(r, r'; \Ei-E f \)\ ^(r') ^/(r), (1) 
where W(r, r'; u) is the frequency-dependent dynamically screened interaction 

W(r, r'; cu) = t>(r, r') + / cfri / <ir 2 f(r, ri) 

xx(ri,r 2 ;^)w(r 2 ,r'), (2) 

v (r, r') and x( r i? r 2;^) being the bare Coulomb interaction and the exact density-response 
function of the Fermi system, respectively. 

In the case of a bounded 3D electron gas that is translationally invariant in the plane of 
the surface, which we assume to be normal to the z axis, the probe-particle initial and final 
states are of the form 

^(r) = i ( 2 )e ik "- r (3) 

and 



M*) = M*)j° H ~ %) ~> (4) 



with energies 



k 2 

Ei = Si + -±- (5) 

2m; 



and 



(k|| -qil) 2 

where the one-particle wave functions <j)(z) and energies e describe motion normal to the 
surface. Using these wave functions and energies in Eq. ([!]), one easily finds 



' = 2E/(|L/d,/d^-( 



Z)0 f {Z 



xlm 



-Wiz^'-^q^lEi-Efl) 



\[z (p f \z 



W(z, z'; qii, ui) being the 2D Fourier transform of the screened interaction. 



3. SCREENED INTERACTION 



(7) 



The main ingredient in the evaluation of electron-hole decay rates in a bounded 3D 
electron gas is the screened interaction W(z, z'; qii, a;). When both z and z' are fixed far 
from the surface, a few atomic layers within the bulk, there is translational invariance in 
the direction normal to the surface, and W(z, z'; a;) can then be easily obtained from the 
knowledge of the dielectric function of a homogeneous electron gas e(q, u): 



(8) 



where v (q) represents the three-dimensional Fourier transform of the bare Coulomb interac- 
tion, and q 2 = q 2 + q 2 - In the random-phase approximation (RPA), e(q, uj) is the Lindhard 
dielectric function ]T7 . 



Shockley surface-state electrons and holes are found to be very close to the Fermi level 
(\E — Ep\ « Ep)- Hence, as the energy transfer \Ei — Ef \ entering Eq. (0) cannot exceed 
the value \E — Ep\, one can use the low-frequency form of the Lindhard dielectric function, 
which yields the energy-loss function 

2 



Im 



-e 1 (q,u) 



,0)uQ(2q F -q), 



(9) 



where Q(x) represents the Heaviside function and qp is the Fermi momentum. If one further 
replaces the static dielectric function e(q, 0) by the Thomas-Fermi (TF) approximation, then 
one finds 



Im 



-e 1 (q,u) 



2q 



(q 2 + q 2 TF ) 2 
5 



ujQ(2q F - q), 



(10) 



Qtf = y^QF/ff being the TF momentum. 



occurs at z = z' and is 



In the bulk, the maximum magnitude of Im —W(z,z';q\\,uj) 
independent of the actual value of z. Fig. 2 shows this quantity for Cu(lll), as obtained 
from Eq. (§) with either the Lindhard dielectric function or the approximated energy-loss 
function of Eq. (|l0|), versus q\\ for fixed values of uj [uj — 0.2, 0.3, 0.4, 0.5 eV]. One clearly 
sees that for the momentum and energy transfers of interest (see Fig. 1), the approximated 
form of Eq. (|10"D yields a screened interaction which is close to that obtained with use of the 
Lindhard dielectric function, small differences being mainly due to the fact that the static 
dielectric function e(q, 0) has been replaced by the TF approximation. Contributions to the 
surface-state decay rate coming from the penetration of the surface-state wave function into 
the solid are, therefore, expected to be well described by Eq. (|10|). 

Nevertheless, coupling of Shockley surface states with the crystal may also occur, either 
through the evanescent tails of bulk states near the surface or within the surface-state band 
itself. These contributions to the surface-state decay rate are both dictated by the screened 
interaction at z points near the surface, where the representation of Eq. (|8|) is not accurate. 

For a realistic description of the screened interaction, we proceed along the lines reported 



m 



lq| . First of all, we compute the single-particle eigenf unctions and eigenvalues of a 



ID model potential This model potential, which reproduces far outside the surface 

the classical image potential, is chosen so as to describe the width and position of the 
energy gap at the F point (ky = 0) and, also, the binding energies of both the Shockley 
surface state at T and the first (n = 1) image-potential induced state. We then introduce 
these eigenfunctions and eigenvalues into Eq. (^), we also use them to compute the non- 
interacting density- response function |20| , and finally solve an integral equation to derive the 
RPA density-response function z'\ qy, uj) and the screened interaction W(z, z'\ qy, a;). 

As the maximum magnitude of Im — W(z, z; qy, a;) /uj near the surface still occurs at 
z ~ z' (only far from the surface into the vacuum, where the Shockley surface-state amplitude 
is negligible, does the maximum of this quantity occur at z ^ z' |2l[]), we choose z — z' 
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and show Im — W(z, z; qn, uj) /urn. Fig. 3, as a function of z going from the bulk to the 
vacuum, for q» = 0.2 Oq 1 (ao is the Bohr radius) and for various values of uj. We find that 
at the surface Im[— VK]/a; is enhanced, for the smallest values of uj, by a factor of 3 relative 
to the bulk. This is mainly due to the strong anisotropy at the surface, which enhances the 
electron-hole pair creation probability j2"2"fl . We also find that for frequencies smaller than 
0.2 — 0.3 eV, the strong anisotropy of the surface does not modify the linear frequency scaling 
characteristic of the bulk screened interaction [see Eq.®]. However, for the larger but still 
small frequencies explored in Fig. 3, the screened interaction at the surface exhibits an 
important enhancement due to corrections to the linear scaling that are not present within 
the bulk. As a result, for the largest values of uj that we have considered, Im[— W]/uj is 
enhanced at the surface by a factor as large as 4. 

Fig. 4 displays Im — W(z, z; qy , uj) /w as a function of z, for uj = 0.2 eV and for various 
values of q\\. While Im[— W]/uj is enhanced at the surface, for the smallest values of q\\, by 
a factor of 3 or 4 in the frequency range uj = 0.2 — 0.5 eV, this enhancement becomes very 
small for the largest values of q« that we have explored. These calculations indicate that 
only transitions with q < 0.4 a^ 1 are affected by the strong anisotropy at the surface. 



4. LIFETIME BROADENING 

We compute the lifetime broadening of Shockley surface states from Eq. (0) with a 
realistic description of the RPA screened interaction, as described in the previous section. 
Although Eq. (]?]) has been derived by assuming that our electron system is translationally 
invariant in the plane of the surface, we account for the potential variation parallel to the 
surface through the introduction of a realistic effective mass into Eqs. (|5|) and (^), and also 
through the introduction into Eq. ([F]) of initial and final wave functions that change along 
the actual dispersion curve of each state. The effective mass of bulk states has been chosen 
to increase from our computed values of 0.31, 0.25, and 0.21 m e at the bottom of the gap 
in Cu(lll), Ag(lll), and Au(lll), respectively, to the free-electron mass m e at the bottom 
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of the valence band. The z-dependent initial and final wave functions <pi(z) and <f)f(z) have 
been recalculated for each value of kn and kn — qy, respectively, as was done at the T point, 
with use of the ID hamiltonian described in the previous section. 

As the relaxation of Shockley surface-state electrons and holes may proceed either 
through the interband coupling with bulk states (3D channel) or through the intraband 
coupling, within the surface state itself, with surface states of different wave vector parallel 
to the surface (2D channel), we also consider the decay of electrons and holes in 3D and 2D 
uniform systems. Introduction of Eq. (|i0|) into Eq. (H) and then Eq. (§) into Eq. (0) with 
both <fii(z) and <f>f(z) replaced by plane waves, yields |23| 



'3D 



8q 2 F 



tan 1 y/nqp 



(E-E F ) 
k 



(11) 



1 + 7rq F _ 

k being the momentum of the excited electron or hole. On the same level of approximation, 
the 2D decay rate is given by f25| , 



'2D 



4tt 



In 



\E-E F \ 
E F 



TF 



1 , 2 <?2D 

H h In 

2 q F 



E — E F 
E P 



(12) 



where qj£ = 2m is the TF screening wave vector in 2D, and m is the electron mass. 



A. Surface-state holes at V 

The partially occupied Shockley surface-state band forms a uniform two-dimensional 
electron gas, with the 2D Fermi energy being given by the band-edge of the surface state, 
i.e., £^=445, 67, and 505 eV for the (111) surfaces of Cu, Ag, and Au, respectively. Our 
full calculation for the intraband linewidth of the band-edge surface-state hole on the (111) 
surfaces of Cu, Ag, and Au is presented in Fig. 5 by full squares, together with the linewidth 
of Eq. ( |i~2"D with E = and the electron mass m chosen to be either the free-electron mass 
(solid line) or the surface-state effective mass (full circles). We find that calculations based on 
a pure 2D electron gas are larger than the actual intraband contribution to the linewidth by 
a factor of ~ 7. This large discrepancy is due to the fact that electron-electron interactions 
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within the actual Shockley surface-state 2D band are strongly screened by the underlying 
3D bulk electron system, thereby reducing the scattering probability 

Separate intraband and interband contributions to the linewidth of Shockley surface-state 
holes at the Y point of the projected bulk band-gap of the (111) surfaces of Cu, Ag, and 
Au are displayed in Table I. The 3D linewidth of bulk free holes with the same energy, as 
obtained from Eq. (|Tl|), is also shown in this table. Differences between our full interband 
calculations and those obtained from Eq. (|ll|) arise from (a) the enhancement of Im[— W] at 
the surface, which increases the linewidth, (b) localization of the surface-state wave function 
in the direction perpendicular to the surface, and (c) the restriction that only bulk states 
with energy lying below the projected band-gap are allowed. Both localization of the surface- 
state wave function and the presence of the band-gap reduce the linewidth, and therefore 
tend to compensate the enhancement of Im[— W] at the surface. In the case of Cu(lll) 
this compensation is nearly complete, thereby yielding an interband linewidth that nearly 
coincides with the 3D linewidth of free holes. However, this is not necessarily the case for 
other materials, depending on the surface band structure. 

The impact of the enhanced Im[— W] at the surface on both interband and intraband 
contributions to the Shockley surface-state hole linewidth at the T point in Cu(lll) is 
illustrated in Fig. 6. In this figure, our full calculations of the interband and intraband 
contributions to the total inelastic linewidth (full circles) are compared to the results we 
have also obtained from Eq. (]?]), but with the actual screened interaction replaced by that 
of Eq. (§) (full triangles). 3D and 2D free-electron gas calculations, as obtained from Eqs. 
(|TT1) and (|T2|), are represented by open squares. We find that the impact of the surface 
anisotropy on Im[— W] (difference between full circles and triangles) is to largely increase 
both interband and intraband contributions to the linewidth. Fig. 6 clearly shows that 
the agreement, in the case of Cu(lll), between our interband linewidth (full circle) and 
that obtained from Eq. ( [TT]) (open square) is due to a fortuitous cancellation between both 
localization of the surface-state wave function and the presence of the band-gap, on the one 
hand, and the enhancement of Im[— W] at the surface, on the other. We also note from 
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this figure that a large contribution from intraband transtions, strongly screened by the 
underlying 3D bulk electron system, is responsible for the large differences between the 3D 
free-electron gas prediction and the experimental results (represented by open triangles), as 
discussed in fll3"] . 

B. Surface-state electrons and holes with ky ^ 

Our full calculation of the total inelastic linewidth of Shockley surface-state electrons and 
holes on Cu(lll), as obtained from Eq. fl7|), is shown in Fig. 7 (solid line with circles) versus 
the surface-state energy Separate interband and intraband contributions to the linewidth 
are also represented in this figure, by solid lines with triangles (interband) and inverted 
triangles (intraband), and the linewidth of bulk free-electrons with the same energy, as 
obtained from Eq. (|TT|), is represented by a dotted line. Measurements of the linewidth of 
excited holes at the surface-state band edge (ky = 0) are represented by open triangles, and 
the STS measurements for surface-state electrons reported by Burgi et al (ky 7^ 0) are 
represented by full squares. 

We find that the fortuitous agreement between the actual interband linewidth and that 
derived from Eq. ([TT|) is present not only for surface-state holes at the V point, but also for 
surface-state electrons and holes near the Fermi level. Hence, the interband linewidth shows 
the (E — Ep)~ 2 energy dependence that is present in the case of a 3D free-electron gas. 

Fig. 7 shows that at the surface-state band edge (ky = the intraband linewidth 
represents an 80% of the total linewidth (see also Table I). Nevertheless, as the surface-sate 
wave vector parallel to the surface increases, the surface-state wave function acquires a bulk- 
like character, with a larger penetration into the bulk, and the intraband contribution to 
the linewidth decreases very rapidly. This conclusion explains the experimental observation 
that while 3D free-electron gas predictions of surface-state hole lifetimes are too large, they 
are in the case of surface-state electrons (above the Fermi level) comparable to measured 
lifetimes. 
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Open triangles in Fig. 7 represent the inelastic contribution to measured linewidths of 
Shockley surface-state holes at the T point, which are to be be compared with our calcu- 
lated linewidth at this point (full circle). Since our model, which correctly reproduces the 
behaviour of the s-p valence states, does not include screening from ci-electrons, differences 
between measured linewidths (open triangles) and our calculations (full circles) are expected 
to be due to the effect of virtual transitions, giving rise to additional screening by (^-electrons 
26fl , as discussed in [ |13|] . 

The STS measurements represented in Fig. 7 by full squares include both inelastic and 
electron-phonon contributions to the surface-state linewidth. Hence, one must be cautious 
in the comparison of these measurements with our calculations. Electron-phonon linewidths 
of typically 8.0 in Cu(lll) |27|, which are essentially independent of the surface-state energy, 
yield theoretical predictions for the total linewidth that are above the experimental obser- 
vation, especially for the lowest surface-state energies explored. This discrepancy between 
theory and experiment could be partially compensated by the screening of (i-electrons, which 
as in the case of surface-state holes at the T point, reduce the inelastic linewidth. 



5. SUMMARY AND CONCLUSIONS 

We have reported theoretical investigations of surface-state electron and hole lifetimes 
in the noble metals Cu, Ag, and Au. We have also presented self-consistent calculations 
of the screened interaction, and have investigated the role that the partially occupied 2D 
surface-state band and the underlying 3D electron gas play in the relaxation mechanism. 

We have reached the conclusion that, while the linewidth of surface-state excited holes 
at the T point of the (111) surfaces of the noble metals is dominated by a 2D decay channel, 
major contributions to the electron-electron interaction of surface-state electrons above the 
Fermi level come from the underlying bulk electrons. This key dependence of the relative 
contribution of intraband transitions to the total decay rate explains the experimental ob- 
servation that while 3D free-electron gas predictions of surface-state hole lifetimes are too 
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large, they are in the case of surface-state electrons (above the Fermi level) comparable to 
measured lifetimes. 
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FIGURES 

FIG. 1. Dispersion of Shockley surface state and bottom of projected band-gap on the (111) 
surfaces of (a) Cu, (b) Ag, and (c) Au. Shaded areas represent areas outside the band-gap, where 
bulk states exist. Relaxation of Shockley surface-state electrons and holes may proceed either 
through interband coupling with bulk states (3D channel) or through intraband coupling, within 
surface state itself, with surface states of different wave vector parallel to surface (2D channel). 

FIG. 2. Imaginary part of bulk screened interaction Im — W bulk (z, z; qii , u>) of Eq. (g), as a 
function of qu for various values of u: 0.2, 0.3, 0.4, and 0.5 eV. Solid and dashed lines repre- 
sent the result of introducing into Eq. (||) either Lindhard dielectric function (solid lines) or the 
approximated energy-loss function of Eq. ([ToD (dashed lines). 



FIG. 3. Imaginary part of the scaled RPA screened interaction Im —W(z,z;q,»,u) /lo of Eq. 
(|7|) , versus z, for q» = 0.2 ag 1 (ao is the Bohr radius) and various values of u: 0.2, 0.3, 0.4, and 
0.5 eV. 

FIG. 4. As in Fig. 3, for uj = 0.2 eV and various values of qn: 0.2, 0.3, 0.4, and 0.5 Oq . 

FIG. 5. Solid line represents linewidth of holes at the bottom (E = 0) of a 2D free-electron 
gas, as obtained from Eq. (^) with m = 1, versus the the 2D Fermi energy. Squares represent 
our full calculation for the intraband contribution to the linewidth of band-edge (ky = 0) Shockley 
surface-state holes on the noble metals Cu, Ag, and Au. 2D free-electron gas predictions of Eq. (12) 
with Ep = 445 meV and m = 0.42 (Cu), Ep = 67meV and m = 0.44 (Ag), and Ep = 505 meV and 
m = 0.28 (Au) are represented by circles. Triangles represent experimentally determined inelastic 
linewidths taken from |0|] (after subtraction of estimated electron-phonon linewidth of 8, 5.2, and 
5.2 meV for Cu(lll), Ag(lll), and Au(lll), respectively. 
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FIG. 6. Band-edge (ky = 0) Shockley surface-state hole linewidths in Cu(lll). Full circles: 
Separate interband and intraband contributions to total inelastic linewidth, as obtained from 
Eq. (|7|) with our full surface calculation of W(z, 2/; q, \Ei — Ej\). Full triangles: Separate in- 
terband and intraband contributions to the total inelastic linewidth, as obtained from Eq. (|7|) with 
W(z, z'; q, \Ei — Ef\) replaced by that of Eq. (|8|). Open squares: 3D and 2D free-electron gas cal- 
culations of interband and intraband linewidths, as obtained from Eqs. (|i"l|) and (12), respectively. 
Experimentally determined inelastic linewidths of [j^] (after extrapolation of ARP linewidth to zero 
defect density) and (jDJ (after subtraction of estimated electron-phonon linewidth of 8 meV) are 
represented by the open triangle and the inverted open triangle, respectively. 

FIG. 7. Scaled inelastic linewidth of Eq. (|7|), t~ 1 /(E — Ep) 2 , of Shockley surface-state elec- 
trons and holes in Cu(lll), as a function of surface-state energy. Total linewidth is represented 
by a solid line with circles. Intraband and interband contributions to the linewidth are repre- 
sented by solid lines with triangles (interband) and inverted triangles (intraband). 3D free-electron 
gas prediction of Eq. (|i~2|) is represented by a dotted line. Experimentally determined inelastic 
linewidths of []|] (after extrapolation of the ARP linewidth to zero defect density) and (af- 
ter subtraction of an estimated electron-phonon linewidth of 8 meV) are represented by the open 
triangle and the inverted open triangle, respectively, as in Fig. 6. Full squares represent STS 
measurements of linewidth (with no subtraction of electron-phonon contribution) of surface-state 
and surface-resonance electrons reported in [14]. 
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TABLES 

TABLE I. Decay rates, in linewidth units (meV), of the Shockley surface-state hole at the F 
point of the noble metals. Decay rates in a 3D free-electron gas of holes with the energy of the 
Shockley surface-state at T are also displayed. 

Surface Energy 3D Inter Intra Total 

Cu(lll) -445 5.9 6 19 25 

Ag(lll) -67 0.18 0.3 2.7 3 

Au(lll) -505 10 8 21 29 
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